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Abstract 

We examine the vibrational stability of three candidate structures for the B24N24 
cage and report their infra-red (IR) and Raman spectra. The candidate structures 
considered are a round cage with octahedral 0symmetry, a cage with S4 symmetry 
that satisfies the isolated square rule, and a cage of Sg symmetry, which combines 
the caps of the (4,4) nanotube, and contains two extra squares and octagons. The 
calculations are performed within density functional theory, at the all electron level, 
with large basis sets, and within the generalized gradient approximation. The verti- 
cal ionization potential (VIP) and static dipole polarizability are also reported. The 
S4 and Ss cages are energetically nearly degenerate and are favored over the 0cage 
which has six extra octagons and squares. The IR and Raman spectra of the three 
clusters show notable differences providing thereby a way to identify and possibly 
synthesize the cages. 

Key words: boron nitride, nanotubes, cages, fullerene, infra-red, ionization 
potential and polarizability 
PACS: 36.40.-C, 73.22.-f, 61.48.-Fc 



* Fax: -hl-202-767-1716 

Email addresses: rzope@alchemy.nrl.navy.mil (Rajendra R. Zope ), 
baruah@dave.nrl.navy.mil (Tunna Baruah ), pederson@dave.nrl.navy.mil 
(Mark R. Pederson), dunlap@nrl.navy.mil (Brett I. Dunlap). 



To be submitted to Chemical Physics Letters 



2 February 2008 



Recently, boron nitride (BN) cages were synthesized and detected by laser des- 
orption time of flight mass spectroscopy [1]. The B24N24 cluster was observed 
in abundance and its structure was proposed to be a round cage structure with 
octahedral (0) symmetry (Cf. Fig. 1), in an analogy to the icosahedral Ceo 
cluster. It consists of alternate BN atoms and is made round by introducing 
eighteen defects: six octagons and twelve squares, in the hexagonal network 
of a single BN sheet. In the icosahedral Ceo fullerene every carbon atom is 
equivalent by symmetry. Similarly every boron or nitrogen atom in this octa- 
hedral structure is equivalent by symmetry. Taking a cue from the similarity 
of Ceo fullerene and the O B24N24 cage it was suggested that the two halves of 
the round octahedral cage can form hemispherical caps for the (4,4) BN nan- 
otubes [2] analogous to Ceo hemispheres capping the (5,5) carbon nanotube 
[3]. Further by Euler theorem, in analogy to the twelve isolated pentagons of 
Ceo, six isolates squares close the alternate BN fuUerenes [4,5]. For the B24N24, 
the resultant cage (Cf. Fig. 1) has S4 symmetry [5]. In a recent study [6] that 
examined seven candidate structures for the B24N24 geometry, the S4 cage 
was found to be energetically favorable over the 0cage in agreement with our 
work. A new, most stable cage structure was found that has Sg symmetry and 
contains eight squares and two octagons in the hexagonal network [6] . The Sg 
cage is related to the B28N28 or B32N32 (4,4) closed nanotubes in Ref. [2] and 
can be derived from B28N28 by eliminating a middle ring of alternating BN 
atoms followed by rotation by 45°. The energy differences between these iso- 
mers, particularly between the Ss and S4 cages are quite small (O.leV or less) 
and thus based on energetics it is not apparent which isomer is most likely to 
be observed in the experiment. Clearly more theoretical calculations or further 
experimental characterization is necessary. For this purpose we have carried 
out density functional [8] calculations on the three candidate structures of 
the B24N24 clusters. We predict the infra-red (IR) and Raman spectra for the 
three isomers which may aid the experimentalists in identifying the structure 
of B24N24 . Along with these spectra we also report the ionization potential, 
and the static dipole polarizability of the three cages. Our calculations are 
carried out at the all electron level by the NRLMOL code [9,10,11], using 
large polarized Gaussian basis sets [12] optimized for these density functional 
calculations, and using the Perdew-Burke-Ernzerhof generalized gradient ap- 
proximation (PBE-GGA) [13]. 

The energetics and the electronic properties of optimized structures are sum- 
marized in Table 1. The calculations predict the Ss cage as the lowest en- 
ergy structure. The S4 isomer is within 0.03 eV while the 0cage is 2.6 eV 
higher. This energy ordering is in accord with earlier hybrid density functional 
(B3LYP/6-31G*) calculations [6]. These energy differences are of same magni- 
tude even with addition of the zero-point energy. The vibrational frequencies 
determined within the harmonic approximation indicate all structures to be 
local minima on the potential energy surface. Further, consistent with earlier 
prediction on the BN cages and nanotubes [14], these cages are characterized 



2 



A E - 2.6 eV 




Fig. 1. (Color online) Two different views of optimized B24N24 cages. The relative 
energies are given with respect to Ss total energy. 

Table 1 

The total electronic energy E (in Hartree), the zero point energy (ZPE), A E (energy 
relative to the energy of Sg cage and excluding ZPE contribution) (in kJ/mol), 
binding energy (BE) per atom (in eV), the vertical ionization potential (VIP) (in 
eV), and the HOMO-LUMO gap (in eV) of the B24N24 cages. The ZPE is half the 
sum of vibrational frequencies. 





E 


ZPE 


AE 


BE 


VIP 


gap 





-1910.96994 


664.49 


2.6 


8.27 


8.4 


4.7 


S4 


-1911.06552 


678.74 


0.03 


8.32 


8.5 


4.9 




-1911.06676 


675.02 


0.0 


8.32 


8.3 


4.6 



by a wide energy gap of about 4.7 eV between the highest occupied molec- 
ular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). 
We note that the HOMO-LUMO gaps obtained within the present PBE-GGA 
generally underestimate the sum of ionization potential and electron affinity. 
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Table 2 

The static dipole polarizability (in A^) of the 0, S4 , and Ss isomers of B24N24 . 
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Table 3 

The IR active frequencies of B24N24 cages. The quantity inside the bracket is the 
intensity (Debye^/amu/A^) the respective absorption. Only IR active modes with 
intensity greater than 3 (Debye^/amu/A^) are given. 



Symmetry Frequency (cm ^ ) 






Ti 


759 (7), 1471 (83) 


S4 


B 


1368 (7), 1423 (47), 1448 (34) 




E 


1353 (7), 1399 (26), 1442 (30) 


Ss 


B 


1432 (87) 




El 


767 (7), 1406 (27), 1461 (45). 



The vertical ionization potential (VIP) is calculated from the the differences in 
total energy of the self-consistent solution of the cluster and its singly charged 
positive ion, at the optimal ionic configuration of the neutral cluster. The 
calculated VIP are also given in Table 1. We note that the VIP reported for 
the 0B24N24 cage in Ref. [2] should be 8.64 eV instead of incorrectly reported 
value of 6.66 eV. The vertical electron affinity calculations can be performed in 
similar fashion. We find that the additional electron is weakly bound (despite 
relatively large basis set used in the calculation its eigenvalue turns out to be 
positive although the total energy is lowered with respect to neutral cluster). 

The static dipole polarizability is an important physical quantity that charac- 
terizes the electric response of the clusters to applied uniform electric field. It 
is calculated from the so called finite field method in which the self-consistent 
problem is solved for the usual Hamiltonian augmented with the field term 
E ■ f. The polarizabihty tensor is built from the self-consistent solutions per- 
formed for various values of electric field E along different directions. The 
nuclei are assumed to be frozen. The calculated values are shown in Table 2. 
The mean polarizability of the three isomers is of similar magnitude with the 
0cage being slightly more polarizable. The structural differences of the iso- 
mers are reflected in the anisotropy of the polarizability. Unlike the perfectly 
round 0cage, the isomers S% and S4 have nonzero anisotropy. 

The calculated IR spectra are presented in Fig. 2. The frequencies and their 
symmetries of IR active modes are listed in Table 3. The IR absorption in- 
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Fig. 2. The infra-red spectra of the three cage structures of the B24N24 . 

tensities are broadened by 6 cm~^ to mimic experimental uncertainties. The 
IR spectra of the three isomers have different structure. The high symmetry 
of the round 0cage results in its IR spectrum with four peaks. It shows a 
strong absorption at 1471 cm~^ and a weak absorption at 759 cm~^ . Other 
peaks around 950 cm~^ and 500 cm~^ have relative intensity less than 2%. 
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Table 4 

The Raman active frequencies of B24N24 cages. (Frequencies with very weak ab- 
sorption are not presented). 



Symmetry Frequency (cm ^ ) 



(7) 


To 


1 




Ai 


428, 780, 790, 901 




E 


131, 910 


S4 


A 


192, 215, 472, 760, 794, 820 




E 


236 




A 


262, 469, 798 




E2 


158, 247 



The strong absorption mode at 1471 cm~^ is three fold degenerate and corre- 
sponds to stretching and a compression of the alternate bond in the octagon. 
An earlier IR calculation [7] on the 0B24N24 cage that used semi-empirical 
modified neglect of diatomic overlap (MNDO) method reports three very dif- 
ferent IR active modes at 1356, 1336, and 772 cm"^ . The IR spectra of S4 
can be broadly classified into two categories of bands in the frequency range 
1353-1448 cm~^ . The first set consists of small peaks around 1350 cm^^ while 
the second set shows three peaks at 1399, 1423, and at 1442-1448 cm~^ . In 
the experiment, the later set of peaks may show up as a broad band beginning 
with a shoulder around 1350 cm~^ . The IR spectra of the Ss structure exhibits 
three conspicuous peaks at 1407, 1432, and 1461 cm~^ . The strongest mode 
at 1432 cm~^ corresponds to a compression of two opposite hexagonal bonds 
accompanied by simultaneous stretching of remaining hexagonal bonds. The 
absorption at 1461 cm~^ is due to the compression and stretching of bonds 
of octagons. A small peak is also observed at lower frequency of 767 cm~^ 
. It corresponds to a vibrational mode that results in puckering of alternate 
atoms. The differences observed in the structure of the predicted IR spectra 
of the three isomers indicate the possibility of identifying the isomer by high 
resolution IR spectroscopy measurement. A sharp, high-frequency peak would 
indicate the presence of the octahedral isomer. 

The Raman spectra of the three B24N24 cages are shown in Fig. 3 and the 
vibrational frequencies at which significant absorption occurs arc presented in 
Table 4. The Raman spectra of three cages show an interesting feature: all 
three cages show an intense absorption in the 420 - 480 cm"^ region. In all 
the three structures this mode is an A mode and corresponds to the breathing 
motion of the cages. The intensity of this peak is comparable in all the three 
structures. The 0cage is marked by an intense double peak in the 131-150 cm^^ 
region which is not seen in the other two spectra. These peaks are associated 
with E and T2 modes. Thus two peaks below 160 cm~^ and another one at 
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Fig. 3. The Raman spectra of the three cage structures of the B24N24 . 



around 430 cm~^ , all with similar intensity will signify a 0structure. On the 
other hand, the S4 structure will show three low Raman peaks in the region 
190 - 230 cm~^ and three others around 800 cm~^ . In experiment, they may 
appear as broad peaks with roughly half the intensity of the most intense 
peak at 430 cm~^ . The 0structure also shows a weak absorption around 800 
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cm^^ . The S8 structure has a low frequency peak around 158 cm~-^ which 
is about half as intense as the peak at 472 cm^^ , and a similar one at 798 
cm~^ . The peak at 158 cm~^ due to a doubly degenerate E mode while the 
one at 798 cm~^ has A symmetry. Further it will show a broad band at 230 
cm~^ which will be missing in the two other structures. These differences in 
the Raman spectra of the 0, S4 , and Ss cages may help in identifying them 
in the experimental spectra. 

To summarize, three candidate structures for the ground state geometry of 
the B24N24 are considered. These are, the octahedral structure proposed by 
the experimentalists, the S4 cage that satisfy the isolated square rule, and 
the Ss symmetric tubule containing two octahedrons and eight squares. The 
harmonic frequency analysis indicate all structures to be vibrationally stable. 
The Ss tubule and S4 are energetically nearly degenerate and are favored over 
the 0cage on the basis of energetics. All the three clusters have wide HOMO- 
LUMO gap and high ionization potential. The infra-red and Raman spectra 
show notable differences and therefore point to possible identification of the 
structure by the IR or Raman spectroscopy. The IR spectra perhaps, as in the 
case of Ceo [15,16], could guide methods for optimizing the production of the 
round cluster. 
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